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\Q We describe a calculation of electromagnetic cascading in radiation and matter in the 

early universe initiated by the decay of massive particles or by some other process. We 
have used a combination of Monte Carlo and numerical techniques which enables us to 

© use exact cross sections, where known, for all the relevant processes. In cascades initiated 

after the epoch of big bang nucleosynthesis 7-rays in the cascades will photodisintegrate 

© 4 He, producing 3 He and deuterium. Using the observed 3 He and deuterium abundances 

we are able to place constraints on the cascade energy deposition as a function of cosmic 
time. In the case of the decay of massive primordial particles, we place limits on the 
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density of massive primordial particles as a function of their mean decay time, and on the 



^ expected intensity of decay neutrinos. 
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1 Introduction 



Electromagnetic cascades in the early universe can be initiated by the decay of massive 
particles [1, 2], or by their annihilation, by cusp radiation of ordinary cosmic strings [3], 
by super-massive particles "evaporating" from superconducting cosmic strings [4, 5, 6], 
by evaporation of primordial black holes, and probably by some other processes. These 
cascades result in the production of 3 He and D by disintegration of 4 He by photons 
in the cascade, and we shall refer to this as "cascade nucleosynthesis". The epoch of 
interest for cascade nucleosynthesis has redshift z <;10 6 , by which time production of 
the light elements D, He and Li by big bang nucleosynthesis has taken place. At this 
epoch, electron-positron pairs are no longer in equilibrium, and black body photons, 7bt>, 
constitute the densest target for electromagnetic cascade development. 

A cascade is initiated by a high energy photon or electron, and develops rapidly in the 
radiation field mainly by photon-photon pair production and inverse Compton scattering: 

e + 7bb^e' + 7', 7 + 7bb^e + + e~. (1) 

Such electron-photon cascading through radiation fields involving photon-photon pair 
production and inverse Compton scattering governs the spectrum of high-energy radiation 
in a variety of astrophysical problems [7, 8, 9, 10, 11, 12, 13, 14, 15]. 

When cascade photons reach energies too low for pair production on the black body 
photons, the cascade development is slowed, and further development occurs in the gas by 
ordinary pair production, but with electrons still losing energy mainly by inverse Compton 
scattering in the black body radiation: 

-f + Z^Z + e+ + e-, e + 7 bb ^e' + 7 / . (2) 

As first pointed out by Lindley [16], since the observed ratios of D/ 4 He and 3 He/ 4 He 
are very small (~ 10~ 5 — 10~ 4 ), cascade nucleosynthesis can put strong constraints on 
the energy going into the particles initiating cascades (7,e + ,e~) in the early universe. 
Cosmological applications of cascade nucleosynthesis have been discussed in several papers 
[16, 17, 1, 2]. The strongest constraints will be placed for redshifts between 10 3 and 10 6 . 
At z > 3 X 10 6 energies of cascade photons are below the threshold for D and 3 He 
production, while at z < 10 3 , direct observation of isotropic X-rays and 7-rays places 
more severe limits on the cascade energy deposition. 

2 The Processes 

Electromagnetic cascades in the early universe take place rapidly in the radiation field, 
and then slowly in the matter. The processes involved in the cascade in the radiation 
field are photon-photon pair production, inverse Compton scattering and photon-photon 
scattering. For interactions in the radiation field, the mean interaction rates for all three 
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processes in black body radiation with temperature T are given by (assuming E ^> mc 2 ) 

r \ f°° r ty m r max ^int(g)(l j i ^ 

T- mt {E,z) = c n[e,T(z)\ d//de, (3) 

where n(e,T) is the differential photon number density of the radiation field, <7; nt (,s) is 
the cross-section for the process, s is the square of the total center of momentum frame 
energy, // = cos 6 is the cosine of the interaction angle, and /3c is the velocity of the 
primary particle (/3 = 1 for photons). Here E is the energy of the primary particle 
and e is the target photon energy. For photon-photon scattering and inverse Compton 
scattering there is no threshold energy, and hence e m i n = and fimax = 1- The threshold 
condition for photon-photon pair production is s > 4m 2 c 4 , giving e m i n = m 2 c 4 /E and 
/"max = 1 — 2m 2 c 4 /(eE). We assume that the present temperature of the microwave 
background radiation is T = 2.735 K. 

Reno and Seckel [18] have explored the consequences of massive particle decay into 
unstable hadrons during the era of primordial nucleosynthesis. Here, particles in the re- 
sulting hadronic cascades interact with nucleons affecting the neutron to proton ratio, and 
hence changing the relative abundances of 4 He, 3 He, D and other light isotopes. For mas- 
sive particle decay at somewhat later epochs (~ 10 3 - 10 7 s), Dimopoulos et al. [19] have 
considered the breakup of 4 He by hadronic cascades. For the epoch under consideration 
in the present paper (> 3 X 10 7 s), unstable hadrons resulting from decay of massive parti- 
cles will decay into neutrinos and an electromagnetic component (electrons and photons) 
before interacting. We therefore only consider the cascade due to the electromagnetic 
component but include in our later discussion the fraction of the massive particle's rest 
mass energy carried away by neutrinos. 

For interactions in the matter, the following processes are included: ordinary (Bethe- 
Heitler) pair production on hydrogen and helium, Compton scattering of energetic photons 
by electrons, photoproduction of pions in photon-proton and photon-helium collisions, and 
bremsstrahlung by energetic electrons on hydrogen and helium. Since the matter density 
depends on epoch as p oc (1 + z)~ 3 } for interactions with matter the interaction rates scale 
with z as 

T^\e,z) = (1 + z) 3 T^\e,z = 0). (4) 

In the case of radiation, the number density of black body photons, rabb = / T)de, 
also depends on epoch as nt,b oc (1 + z)~ 3 . However, because the target photon energies 
also depend on epoch, e oc (1 + z), for interactions with radiation, 

r£?(£, *) = (! + -) 3 lt?{(l + z)E, z = 0}. (5) 

Mean interaction rates are illustrated in Fig. 1 for all of the processes discussed above. 
In the case of interactions with radiation, these are given for the present epoch, (1 + z) = 1, 
but they may be scaled to any epoch shifting the corresponding curve by a factor (1 + z) 
towards lower energies as described in Eq. (5). 
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3 Qualitative Treatment of Electromagnetic Cas- 
cades at Large Redshift 

The black body radiation at the epochs of interest (z < 10 7 ) is characterized by temper- 
ature, T(z) = 2.735(1 + z) K, and photon number density, 

n hh = 4.22 x 10 2 T 2 3 75 (1 + zf cm" 3 , (6) 

where T2.75 = 0.995 is the temperature in units of 2.75 K. 
The density of baryonic gas is 

p b (z) = 1.88 x 10" 29 (1 + zftt b h 2 gem" 3 (7) 

where = Pb/ p c is the mass fraction of baryons in the universe and h is Hubble's constant 
in units of 100 km s -1 Mpc -1 . From the recent review by Copi, Schramm and Turner 
[20], for standard nucleosynthesis one has 

0.009 < tt b h 2 < 0.02, (8) 

and 

0.4 < h < 1.0. (9) 

The baryonic gas consists mainly of hydrogen (~ 77% by mass) and helium (~ 23% by 
mass). The radiation length for gas of this composition is X 66.6 g cm -2 . 

The characteristic interaction rates for photon-photon pair production (PP) and or- 
dinary pair production, i.e. Bethe-Heitler (BH) process, are 



r PP (£, z) « 2.2 x 10" 12 (1 + z f^^-expi-E^/E) s"\ (10) 
at E <^ E a (z) } where 

, , m e c 2 1.12 x 10 6 ^ 

a[z) = ~§F = ~(TT7j~ GeV ' 

and, using the asymptotic pair production cross sections, 

r BH (£,^) w 1.7 x 10" 22 (1 + zf s _1 . (12) 

We must compare these with the expansion rate 

f 3.24 x 10- ls h(l + zf' 2 s" 1 for K2.3x WVth 2 T^ 
[Z) ~ \ 1.65 x 10- 2O T 2 : 2 5 (l + z)- 2 s" 1 for o2.3x iowt 2 : 7 4 5 [ i] 

Strictly, we should use the correct energy-dependent cross sections (described later), but 
this is sufficient for the present discussion. At high energies, the cascade develops entirely 
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on the black body photons by photon-photon pair production and inverse Compton scat- 
tering. The characteristic interaction rates for these processes are much higher than the 
expansion rate, H(z), and thus one can assume the cascade spectrum is formed instantly. 
We shall refer to this spectrum as "the zero-generation spectrum". 

The zero-generation spectrum extends up to a maximum energy, Ec } which is deter- 
mined at low redshifts (z < 10 3 ) by 

T PP (E,z)^H(z), (14) 

and at high redshifts (z > 10 3 ) by 

r P p(£,z)«r B H(£,z). (15) 

The maximum energy obtained in this way is given approximately by 

, , / 4.5 x 10 4 (1 + z)- 1 GeV for z < 10 3 , , 

Ec{z) ~ | 4 ? x 1Q4(1 + z yi Ge y for z > 1Q 3_ (lb) 

For small z, the zero-generation spectrum for a cascade of primary energy E can 
be approximately calculated analytically as described by Berezinsky et al. [11] with the 
result 

f KoiE/Ex)- 1 - 5 for E < E x 
n^(E) « I K (E/E X )- 2 for E x < E < E c (17) 
{ for E c < E 

Here n^(E) gives the differential photon spectrum, Ex = 1.78 X 10 3 (1 + z) 1 GeV is the 
power-law break-energy appropriate for a black body target photon spectrum, and 

En (18) 



u E\[2 + HE c lEx)] 
is the normalization constant. This spectrum is confirmed by Monte Carlo simulation at 

2 = 0. 

When the cut-off energy, Ec given by Eq. (16), is less than the threshold for pho- 
todisintegration of 4 He nuclei (~ 20 MeV), cascade nucleosynthesis is inefficient. This 
condition restricts the epoch of cascade nucleosynthesis to z <Z2 X 10 6 . At large redshifts 
the spectrum is considerably distorted by 77 — > scattering, as was first demonstrated 
by Svensson and Zdziarski [21], and we take this into account in our accurate calculations. 

At redshifts between 10 3 and 2 X 10 6 , where the subsequent cascade development is 
via ordinary pair production and inverse Compton scattering, the zero-generation photons 
survive for a time determined by either the energy loss rate for Compton scattering or 
the interaction rate for ordinary pair production in the gas. During this time they can 
produce light nuclei by photodisintegration. The electrons and positrons give rise to first 
generation photons as a result of inverse Compton scattering. These first generation 
photons then produce the second generation photons, etc. Each generation of photons is 
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strongly shifted to low energies because the inverse Compton scattering is in the Thomson 
regime, and only 1-2 generations of photons are sufficiently energetic to induce cascade 
nucleosynthesis . 

At redshifts z <;10 3 , when interaction times become larger than the Hubble time, only 
the zero-generation photons are produced, and the effectiveness of cascade nucleosynthesis 
diminishes as z decreases. At these redshifts, however, cascade photons can be observed 
directly and such direct observations more strongly constrain the energy density of massive 
primordial particles. 



3.1 Approximate calculation of cascade nucleosynthesis 

For the redshift range 10 3 - 2 X 10 6 the zero generation photons survive ordinary pair 
production and Compton scattering for a time [Tbh^) + k(E)Tcs(E)]~ 1 where 

k ^* l - M2E/3ni/2 (19) 

is the average fraction of energy lost in Compton scattering. Neglecting subsequent gen- 
erations, the number of D nuclei produced by the cascade is given approximately by 

N D « n He c J [T BH (E) + k(E)T cs (E)}- 1 n^(E)a^(E)dE (20) 

where «H e is the number density of He nuclei, <td is the effective cross section for photo- 
disintegration of He into D. 

The effective cross section is the sum over partial cross sections of all channels giving 
rise to the nucleus in question, weighted by the multiplicity. For example, for photodis- 
integration of 4 He into D we have 

<rf (£) = ff ( 7 ,pnD;£) + 2 ff ( 7 ,DD;£), (21) 

and for photodisintegration of 4 He into 3 He we have 

<r% e (E) = (7(7, 3 He n- E) + <r( 7 , 3 H p; E) (22) 

where we have included production of 3 H because it decays into 3 He. For the photodis- 
integration cross sections, we have used data of Arkatov et al. [22]. The effective cross 
sections used for photodisintegration are plotted in Fig. 2(a). We see that the important 
photon energy range is between 25 MeV and 100 MeV and that, above threshold, the 
cross section for production of 3 He is much higher than for production of D. 

In the energy range where photodisintegration of 4 He is important, the pair production 
cross sections have not yet reached their asymptotic values, and are strongly energy 
dependent. We show in Fig. 3 the pair production cross sections for hydrogen and helium 
for the two cases: fully ionized matter, and neutral matter. For E < 100 MeV, the 
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cross sections are almost independent of ionization state, but are quite different from the 
asymptotic values. The interaction rate for pair production is given by 



Fbu(E) = [n H cr lH (E) + nu e cr lH e(E)]c. 



(23) 



Hence, from Eq. (20) we obtain 



I 



(Y/4)og(E)n<?)(E) 



dE 



(24) 



(f - Y)a lH (E) + (Y/4)a lHe (E) + (1 - Y/2)k(E)a C s{E) 



where Y is the fraction of helium in the early universe by mass, and acs(E) is the Klein- 
Nishina cross section. A similar equation governs the number of 3 He nuclei produced. 

4 Accurate Calculation 

To take account of the exact energy dependences of all the cross sections we use the 
Monte Carlo method. However, direct application of Monte Carlo techniques to cascades 
dominated by the physical processes described above over cosmological time intervals 
presents some difficulties, which we will try to address in the following sections. 

The approach we use here is based on the matrix multiplication method described 
by Protheroe [9] and subsequently developed by Protheroe & Stanev [14]. We use a 
Monte Carlo program to calculate the yields of secondary particles due to interactions 
with the thermal radiation and matter. The yields are then used to build up transfer 
matrices which describe the change in the spectra of particles produced after propagating 
through the radiation/matter environment for a time St. Manipulation of the transfer 
matrices as described below enables one to calculate the spectra of particles resulting 
from propagation over arbitrarily large times. 

4.1 Matrix method 

We use 110 fixed logarithmic energy bins of width A log E = 0.1 covering the energy range 
from 10~ 3 GeV to 10 8 GeV. For example, the energy range of the jth energy bin runs from 
]_OC?' -31 )/ 10 GeV to lO^ -30 )/ 10 GeV. The energy spectra of electrons and photons in the 
cascade at time t are represented by by vectors Fj{t) and Fj(t) which give respectively 
the total number of electrons, and photons, in the jth energy bin at time t. 

The numbers of nuclei produced by photodisintegration of 4 He nuclei by photons in 
the cascade are also represented by vectors, Fjit) and Fj(t), which give respectively the 
total number of 3 He nuclei and D ( 2 H) nuclei produced by interactions of photons having 
energy in the jth energy bin at time t. That is, in this case the energy bin index refers to 
the energy of the photon responsible for the photodisintegration, and not to the energy 
of the produced nucleus, which is negligible. 

We define transfer matrices, Tf^(St) } which give the number of particles of type v = e 
(electron), 7 (photon), 3 ( 3 He) or 2 (deuterium) in the bin j which result at a time St 



7 



after a particle of type // =e or 7 and energy in the bin i initiates a cascade. Then, given 
the spectra of particles at time t we can obtain the spectra at time it + St) 



110 



FJ(t + St) 



F-'(t + St) 



Ff(t + St) 



Ff(t + St) 



i=j 
110 

i=j 
110 

i=j 
110 

E[ T ?(W ? (*) + ^(^(*) 



(25) 
(26) 
(27) 
(28) 



where F?(t) and F^it) are the input electron and photon spectra (number of electrons or 
photons in the zth energy bin). 
We could also write this as 



[F(t + St)] = [T(St)][F(t)] 



where 



[F] 



4.2 Transfer matrix calculation 
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(29) 



(30) 



The transfer matrices depend on particle yields, Y"^\ which we define as the probability 
of producing a particle of type f3 in the energy bin j when a primary particle with energy 
in bin i undergoes an interaction of type a. To calculate we use a Monte Carlo 
simulation. For inverse Compton scattering and photon-photon pair production we have 
used the computer code described by Protheroe [9, 12], updated to model interactions with 
a thermal photon distribution of arbitrary temperature. For photon-photon scattering, 
we have used the cross sections given by Berestetskii et al. [23]. 

In the case of inverse Compton scattering in the Thomson regime, the basic matrix 
method fails to predict correctly the electron spectrum, and hence the emitted photon 
spectrum. This is because the fraction of energy lost per interction is small, and effectively 
the electrons suffer continuous energy losses, 



dE 
dt 



-bE 2 



(31) 



Injection of one electron with energy E at t = should result in one electron with energy 

E(t) w (bt + E- 1 )- 1 (32) 
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at time t (i.e. there is very little spread in the final energy). However, the basic matrix 
method would give rise to a broad energy distribution with mean energy equal to Eit). 
In the present problem, electrons in the Thomson regime lose energy by inverse Compton 
scattering at a rate very much higher than by the competing process, bremsstrahlung. 
Therefore we immediatly replace any electron produced with energy in the Thomson 
regime, or in the transition region between Thomson and Klein-Nishina regimes, by all 
the photons from inverse Compton scattering that would be produced while the electron 
subsequently cools. We define a matrix Gf^ which gives the number of photons produced 
in energy bin j by an electron injected with energy in bin i < m. We have chosen the 
maximum energy of electrons treated in this way to be just below the photon-photon pair 
production threshold, i.e. m = 77 at the present epoch. For 61 < i < 77 (at the present 
epoch) the Monte Carlo method is used to calculate Gf^ for j ' > 61, and the remaining 
array elements, for which the electron is well inside the Thomson regime, are obtained 
from the distribution function which we obtain by numerical integration, 

dN = 9m e c 2 Jo 7 J77- 4 / &£- 1 n(e)/(£ 1 ) 

dE 16fdeen(e) ' 1 ' 

where /(£i) is given in ref. [24]. 

For interactions with matter (ordinary pair production, bremsstrahlung, Compton 
scattering, photoproduction) we assume a composition 77% hydrogen and 23% helium 
by mass. Photon and electron yields from photoproduction follow from the 7r ± — — > 

and ir° — > 27 decays. In the case of ordinary pair production and bremsstrahlung, 
yields are calculated for both neutral and fully ionized matter, and mixed in a ratio 
appropriate to the fractional ionization at the epoch for which the transfer matrices are 
to be calculated. For the fractional ionization we use Eq. (3.95) of Kolb and Turner [25]. 

It is important to emphasize that in the most important energy range (20 - 100 MeV) 
the cross-sections for bremsstrahlung and pair production are very strongly energy de- 
pendent. For pair production we use the direct calculations of Hubbel, Gimm and Overbo 
[26] for hydrogen and helium, in which the screening correction (for neutral matter) is 
explicit. For bremsstrahlung we use the expressions of Koch and Motz [27] with form 
factors for hydrogen and helium adjusted to represent the more precise values of Tsai 
[28]. The cross-sections for fully ionized matter are calculated with the bremsstrahlung 
formulae valid in the absence of screening by the atomic electrons. 

In the case of production of 3 He and deuterium nuclei, the yield is simply the effective 
cross section for production of the nucleus in question divided by the total cross section. 
There is, however, a further complication. During pion photoproduction on 4 He the nu- 
cleus almost always fragments, and we must therefore take account of photodisintegration 
during pion photoproduction. We know of no experimental data giving branching ratios 
for the various possible final state nuclei, and new measurements are urgently required. 
In the mean time, we assume this process to be similar to the break-up of 4 He during 
collisions with nucleons in which pions are produced. Following Meyer [29] we assume 
final states ( 3 Hen) : ( 3 Hp) : (DD) : (pnD) are produced in the ratio 2:2:1:2 with neg- 
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ligible production of the final state (2p2n). Hence the effective cross sections used for 
photodisintegration during pion photoproduction are 

4 

cr e ff (D, pions) = cr e ff( He,pions) = — cr(photoproduction). (34) 

The effective cross sections used for photodisintegration during pion photoproduction 
have been added to Fig. 2(a). In Fig. 2(b) we have plotted the effective cross sections 
for both photodisintegration and disintegration during photoproduction multiplied by 
E~ ' 7 . For this log-linear plot, the areas then show the relative contributions to 3 He 
and D production for an E~ x ' 7 photon spectrum (a single power-law approximation to 
the photon spectrum in the cascade). Clearly, disintegration during photoproduction is 
particularly important for calculating abundances of deuterium, but is less important for 
3 He production. Since the inferred ratio of ( 3 He + D) to 4 He is of the order of 10~ 4 , 
this implies that photodisintegration of 3 He and D is unimportant in determining their 
abundances. We therefore neglect this process. 

To calculate the transfer matrices we have used a modification of the semi-analytical 
technique described by Protheroe & Stanev [14]. From Fig. 1 we see that at all epochs 
the highest interaction rate is that for inverse Compton scattering by electrons at low 
energies where the scattering is in the Thomson regime, Tic — > Tj = n^crxc. If St is 
much shorter than the shortest interaction time in the cascade, i.e. St <C 1/Tt, then 

T^ist) « %[i-^r e (^0] + ^[ric(^0^ ce + r br e m (^)^ reme ], (35) 
T:](st) « a[r IC (E,-)^ C7 + r brem (E,-)^ rem7 ], (36) 

T]USt) « 8t[T FF (Ei)Y? p e + T cs (E t )Y l f e + T BU (E t )Y l B l s e + T photo (E t )Y? hot ° e ] } (37) 



T!]{St) « S^l-StY^Ei)}^ 



+St[VUEi)Y:p tl + T cs (E t )Y^ s ~< + r photo (£0^-f otO7 ]> (38) 
T<*(St) « 0, (39) 



Tf 3 (St) « 0,, (40) 

T^(st) « ^[r P 4iw 3 +n*±(^r o3 ]> (4i) 

Tg(6t) « StiV^E^^ + V^(E i )Y^°\ (42) 
T^(St) = (43) 

T?/(ft) = % (44) 
where 

T e (Ei) = T IC (Ei) + r brem (£ 8 ), (45) 

and 

r 7 (E,-) = r PP (E t -) + r 8cat (E,-) + r cs (E,-) + r photo (E,-) + r BU (E t ) + v PD (E t ) (46) 

give the total interaction rates of electrons and photons. In the equations above, we 
have used the following abbreviations: IC (inverse Compton), brem (bremsstrahlung), 
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PP (photon-photon pair production), CS (Compton scattering), BH (Bethe-Heitler pair 
production), scat (photon-photon scattering), photo (photoproduction on both hydrogen 
and helium), photo (He) (photoproduction on helium), and PD (photodisintegration). 

Modifications outlined earlier to take account of inverse Compton scattering properly 
in the Thomson regime and the transition region between the Thomson and Klein-Nishina 
regimes (i.e. for j < m) are as follows: 

m 

vim - vjm + Y.mcimir + ^uE^Yt-^, (4?) 

k=i 

m 

T%(6t) <- TZ(6t) + ^6t\r„(E i )Y? k *° + T ca (E i )Y£ s ° + 

k=i 

+r BU (E i )Y™ e + r photo (i^r o e ]cg (48) 

T^(«) = SiAl - StT e (Ei)] (49) 
T£(6t) = 0. (50) 

We require I /St be much larger than the largest interaction rate in the problem, 
Tt ~ 10 — 11 (1 + z) s -1 , and hence typically we use St 10 lo /(l + z) s. The cascade is 
followed for a time t max which must be much longer than the largest interaction time for 
interactions with matter at the energies at which photodisintegration and disintegration 
during photoproduction can be significant, say from 30 MeV to 10 GeV. At about 10 GeV, 
Tcs ~ 2 X 10~ 25 (1 + z) s -1 so we would require t max ~ 5 X 10 25 /(1 + z) s. To complete 
the calculation of the cascade over time t max using repeated application of the transfer 
matrices would therefore require t meix /6t ~ 5 X 10 15 steps. This is clearly impractical, and 
we must use the more sophisticated approach described below. 



4.3 Matrix doubling method 

The matrix method and matrix doubling technique have been used for many years in 
radiative transfer problems [30, 31]. The method used here to calculate the spectrum of 
particles emerging after an arbitrary time is that described by Protheroe & Stanev [14], 
and is summarized below. Once the transfer matrices have been calculated for a time St, 
the transfer matrix for a time 2St is simply given by applying the transfer matrices twice, 

[T(2St)] = [T(St)r. (51) 

In practice, it is necessary to use double precision and to ensure that energy conservation 
is preserved after each doubling. The new matrices may then be used to calculate the 
transfer matrices for a time interval ASt, and so on. A time interval 2 n St only requires 
the application of this 'matrix doubling' n times. The spectrum of electrons and photons 
after a large time interval At is then given by 

[F(* + A*)] = [T(A*)][F(*)] (52) 
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where [F(t)] represents the input spectra, and At = 2 n 6t. In this way, cascades over 
long time intervals can be modelled quickly and efficiently. For example, to simulate the 
cascade over t max ~ 5 X 10 25 /(1 + z) s with an initial step size of St ~ 10 lo /(l + z) s would 
take only 52 steps. 

As a test we have run the program over such large time intervals and switched off all 
processes except photon-photon pair production, inverse Compton scattering and photon- 
photon scattering so that our results could be compared directly with those of Svensson 
and Zdziarski [21]. In their calculations, Svensson and Zdziarski continuously inject into a 
radiation field photons or electrons with energies above the threshold for photon-photon 
pair production, and solve the kinetic equation to find the steady-state spectrum, for 
some given constant escape time t esc , which they refer to as the "escaping spectrum". 
Two important energies enter in the problem, the maximum photon energy E m defined 
as the energy at which the photon-photon pair production rate equals the escape rate, 
and energy E c defined as the energy at which the photon-photon scattering rate equals 
the escape rate: 

T PP (E m )=t£, (53) 

r scat (£ c ) = C. (54) 

For the case of a black body radaition field Svensson and Zdziarski show results of their cal- 
culations, both with and without photon-photon scattering, for the case where E m / E c = 4. 
In order for us to compare with these results, it was necessary first to work out the value 
of t esc such that E m / E c = 4 for black body radiation of a given temerature. The matrix 
program was then run for injection of a primary electron of energy E{ ^> E m to obtain the 
spectum of photons time t after injection, F 1 (E } t) } for < t < 10t esc . We then obtain 
the escaping spectrum per energy injected by integration, 

1 fimu / f 

F IAE) = V — Fi(E,t)exp(- — 

-^W'esc ^ ''esc 

In Fig. 4 we plot r)(E) = E 2 FJ SC (E) against E / E m for our calculations both including and 
neglecting photon-photon scattering. Our results are compared with those of Svensson and 
Zdziarski [21]. Note that the results of Svensson and Zdziarski are arbitrarily normalized 
such that the no-scattering case has rj(E m ) = 1. We see that the agreement is satisfactory. 
We also show the spectrum (Eq. 17) used in our approximate treatment and the spectrum 
given in Equation 11 of Ellis et al. [1]. 

5 Redshifting 

For (1 + z) <;10 5 , we must take account of redshifting because the expansion rate of the 
universe, H } becomes comparable to the interaction rates in matter. The approach we 
adopt is to propagate initially over a time interval At which would give rise to a change 
in log(l + z) equal to the width of the energy bins, 

Alog(l + z) = 0.1. (56) 




12 



We then redshift the energy bin contents of the vectors representing the photon and 
electron spectra: 

F?<-F? +1 ; F?<-F? +1 . (57) 

Much of the cascade development takes place in this first interval. A further application 
of the transfer matrices for further propagation over an additional time At would give 
rise to a change in log(l + z) approximately equal to the width of the energy bins. We 
then redshift the energy bin contents of the vectors representing the electron and photon 
spectra as described above. While not strictly exact, because the further redshift change 
does not correspond exactly to Alog(l + z) = 0.1, we have found that the error induced by 
this procedure is insignificant in the present problem. This procedure is repeated until the 
cascade is complete, and because the redshifting means that particle energies are more 
rapidly reduced below interaction thresholds, the cascade finishes earlier than without 
redshifting. 

6 Results and Discussion 

We have performed the cascade calculation to find the number of deuterium and 3 He 
nuclei produced by a cascade initiated at the epoch of redshift z c . We give the number of 
nuclei per 1 GeV of total cascade energy, so that the total number of nuclei is obtained by 
multiplying by the total energy of the cascade in GeV. Because of the almost instant for- 
mation of the zero-generation spectrum, the exact shape of the 7-ray or electron injection 
spectrum is of no consequence for further cascade development, and only the total amount 
of injected energy is relevant. The results are given for redshifts at which cascades are 
initiated, z c , in the range 10 2 - 10 7 , and for the following values of h (and 0.4 (0.125); 
0.7 (0.025); and 1.0 (0.01). The resulting number of 3 He and D produced per GeV of the 
total cascade energy is given in Fig. 5(a) and Fig. 5(b) respectively. At redshifts less than 
~ 10 5 for which the maximum energy of photons in the cascade exceeds the thresholds for 
production of 3 He and D by photodisintegration of 4 He, approximately ten times as many 
3 He nuclei are produced compared with D nuclei. This is to be expected given the cross 
sections for photodisintegration and disintegration during photoproduction (see Fig. 2). 

In Fig. 6 we show the sensitivity of the results to: assumptions about the effective 
cross sections for photodisintegration; use (or neglect) of the energy energy dependence 
of the pair production and bremsstrahlung cross sections; and inclusion (or neglect) of 
photon-photon scattering. In this comparison, we use the h = 0.7 results. The full 
curve is the accurate calculation including all effects. The dotted curve shows the effect 
of neglecting photon-photon scattering, which is seen to have the greatest effect on the 
results for z c > 3 X 10 4 . The dashed curve shows the effect of neglecting the disintegration 
of 4 He during pion photoproduction and is seen to be negligible for 3 He production, but 
accounts for up to 50% of all deuterium production, depending on z c . The dot-dash curve 
shows the effect of neglecting the disintegration of 4 He during pion photoproduction, using 
asymptotic pair production and bremsstrahlung cross sections, and neglecting photon- 
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photon scattering. Finally, The heavy dot-dot-dot-dash curve gives the result of the 
approximate treatment. Here, we simply apply Eq. (24) using the energy-dependent cross 
sections, include disintegration during photoproduction, and use the photon spectrum 
given by Eq. (17). 

We are now in a position to summarize the cascade nucleosynthesis scenario in the 
context of the results given in Figs. 5 and 6. The epoch of cascade nucleosynthesis is 
limited by z max and z m i n . The maximum redshift is determined by the condition that the 
maximum energy of photons in the cascade spectrum, given by Eq. (17), must be larger 
than the threshold for D or 3 He production on 4 He, E t h ~ 20 MeV. This condition results 
in 



which is clearly observed in Figs. 5 and 6. 

These Figures demonstrate also that the effectiveness of D and 3 He production de- 
creases as z decreases. The reason is not the decrease in the density of 4 He, as one naively 
might think, but rather the decrease in the number of low-energy photons in the cascade. 
In fact, the lower the redshift, the higher the photon energies in the cascade (both Ex 
and Ec increase), and therefore a smaller fraction of the photons participates in cascade 
nucleosynthesis. One can understand this in a semi-qualitative way from Eq. (20). If we 
neglect the energy dependence of Tbh, and use n^\E) given by Eq. (17), we obtain 



In deriving this Equation, we took into account that, at small z, only the low-energy part 
of the spectrum (Eq. 17) effectively takes part in nucleosynthesis (Ex is high). From 
Eq. (59) we see that the factor (1 + z) 3 in riYi e (z) is compensated for by the same factor in 
Pb(z) } and then we are left with a factor (1 + z)~ ' 5 coming from Ex(z). Numerically, the 
(1 + z)~ ' 5 dependence for the curves in Figs. 5 and 6 would give rather a bad fit because 
it neglects the dependence of Tbh on E and the cascading in the gas. 

At z < 10 3 the cascade photons can be directly observed, and the upper limit for the 
isotropic gamma-ray flux at 10 — 200 MeV is more restrictive for the cascade production 
than nucleosynthesis. Therefore, the most effective epoch for cascade nucleosynthesis 
corresponds to redshifts in the range 10 3 — 2 X 10 6 . 

From Fig. 6 one can see that that role of 77 — > scattering is important only for 
epochs with redshifts z > 5 X 10 4 . This is easy to understand from Fig. 1 which shows 
that, at z = 0, r scat is practically everywhere below the curves ^p — > ee (ordinary pair 
production) and 77 — > ee (photon-photon pair production). For the epoch with redshift 
z we must coherently shift both curves (77 — > and 77 — > ee) by a factor (1 + z) to the 
left. At redshifts z > 10 5 the crossing point of the 77 — > and ^p — > ee curves is at E < 
1 GeV where the ^p — > ee interaction rate is less because of the energy dependence of 
the cross section near threshold, and so 77 scattering becomes relatively more important. 
Passing through the energies where 77 — > scattering dominates, the spectrum of 
the cascade changes, the main effect being that the cut-off energy in the spectrum is 



max 



2.4 x 10 b , 





(59) 
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lowered. The scattered photons do not interact again with the target photons; they are 
just redistributed over the spectrum producing a small bump before the cut off. Both 
these effects (bump and early cut-off) are prominent only at high redshifts z > fO 5 for 
the reasons given above. 

We distinguish the zero-generation cascade from the cascades of the first, second, etc., 
generations. The zero-generation cascade develops on the black body photons mainly 
due to 7 + 7bb — > e + + e~ and e + 7bb — *■ e + 7 scattering. The characteristic times for 
these processes are much shorter than for all other processes, and one can assume that 
the zero-generation spectrum is formed instantly. At small z the spectrum is given by 
Eq. (17); at large z it is distorted by 77 — > scattering, and by energy dependence of 
the inverse Compton scattering. 

At large redshift, z ~ £ max , only photons from the zero-generation cascade participate 
in cascade nucleosynthesis. In this case, the maximum cascade energy, Ec(z) } is close to 
the threshold of the nuclear reactions. The photons of the first generation are strongly 
shifted towards low energies, and thus they become sterile. From Fig. 6 one can see that 
the approximate calculation (dot-dot-dot-dash curve), with the zero-generation photons 
only, is extremely close to the exact calculation in which 77 — > scattering is neglected 
(dotted curve). 

The cascade in the gas develops due to 7 + Z — > Z + e + + e~ and e + 7bb — *■ e + 7 
scattering. In the latter process, the scattered photon is strongly shifted to lower energies 
in comparison with the initial electron because the inverse Compton scattering is in the 
Thomson regime. As a result, at z ~ 10 3 only two generations of photons are sufficiently 
energetic to induce cascade nucleosynthesis. On the other hand, for cascading in the gas 
we require 

J ° p b cdt > X (60) 

where X = |X is the interaction length in g cm -2 (X is the radiation length). This 
implies 

(3H X \" 3 , / 0.01 h \ 2 ' 3 

^(wa) =1 ' 2x1 ° (sw?^J • (61) 

Therefore for z <C 10 3 , again only the zero-generation photons take part in nucleosynthe- 
sis. 



7 Applications and Conclusions 

A few words about some applications of our results are now in order. Cascade nucleosyn- 
thesis strongly restricts high energy processes at 10 3 < z < 10 6 . From Fig. 5, the number 
of 3 He and D nuclei produced by each cascade with total energy E (in GeV) is A^( 3 He) ~ 
(0.1 - l)E and N(B) ~ (1 - 6) x lO" 3 ^, where E is in GeV. The total production of 
3 He and D during the cascade nucleosynthesis must be less than the observed primordial 
quantities. 
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The physical processes ol interest include the decay of heavy relic particles [1, 2], the 
cusp radiation of cosmic strings [3], and the massive particle production by superconduct- 
ing cosmic strings [4, 5, 6]. 

We shall give here two examples of applications. In the first example, we obtain the 
limit on the density of long-lived particles X which can decay into a cascade-producing 
particle c (X — > c+ anything) with a branching ratio b c . It is easy to calculate the 
fraction ( 3 He + D)/H at z = as 

3 He + D /A!! x/^\ [jV(3HeiTx) + iv(DiTx)] ((j2) 



H J n H \lGeV / 

where n* x is the space density that X particles would have at z = if they were stable, 
all densities njj and n* x are taken at z = 0, f c is the fraction of mass mx transferred to 
cascade energy, Tx is the lifetime of the X particle and 

iV( 3 He, T X ) = — [ iV( 3 He, z c (t j) exp(-t/T X )dt. (63) 

Here iV( 3 He, z c ) is the number of 3 He nuclei produced per GeV of the total cascade 
energy, given in Fig. 5(a). iV( 3 He,Tx) obtained by convolving the data of Figure 5 with 
an exponential distribution of decay times with mean decay time Tx (Eq. 63) is given in 
Fig. 7(a). iV(D,Tx), given in Fig. 7(b), is obtained from N(D,t(z c )), given in Fig. 5(b), 
in exactly the same way. 
Defining 

^x = - JLJL , (64) 
Pc 

assuming baryonic matter is 77% hydrogen by mass, and rearranging Eq. (62) we obtain 

In Fig. 8 we use the results of Fig. 7 for the total 3 He plus D produced per GeV of cascade 
energy, together with the upper limit inferred from measurements of 3 He in meteorites 
and the solar wind making assumptions about stellar processing and galactic chemical 
evolution[20] ( 3 He + D)/H < 1.1 X 10~ 4 , to obtain an upper limit to Q* x (tx). Note, 
however, that very recent measurements of D/H are closer to 2.5 X 10~ 4 [32, 33], and if 
such higher values for ( 3 He + D)/H were adopted, the upper limits we would derive would 
be correspondingly higher. Since O^(rx) is proportional to ftb/f c b C} we have plotted 
f c b c ft* x (Tx)/flb- In Fig. 8 we also plot the result of Ellis et al. [1] and a result we 
would have obtained if we had uses the asymptotic Bethe-Heitler pair production and 
bremsstrahlung cross sections. The discrepancy between our results and those of Ellis et 
al. can be partly explained in terms of the cross section. Note that O^(rx) depends on 
( 3 He + D), i.e. mainly on 3 He which is insensitive to assumptions about 4 He disintegration 
during pion photoproduction. From Fig. 8 we find that for mean decay times Tx ranging 
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from 1 year to 10 6 years, a density of dark matter in the form of massive primordial 
particles of only 0.1% to 0.3% of that of normal matter could account for all of the 
observed 3 He and deuterium. 

As a second example, we obtain an upper limit for the neutrino intensity produced by 
the decay of X — > v + anything with a branching ratio b v . We shall assume, as in the 
first example, that these particles decay also into cascade-producing particles (X — > c + 
anything) with a branching ratio b c . Since each X particle decay results in b v neutrinos, 
the present density (z = 0) of neutrinos is n v = b v n* x . Putting n* x = b~ 1 n l/ into Eq. (62) 
and using ( 3 He + D)/H < 1.1 X 10~ 4 one obtains an upper limit for the integral neutrino 
intensity I v = (c/47r)n J/ , 

h < 2.3fM 2 -^- (^^) [iV( 3 He, Tx ) + N(D, t x )]-\ cm" 2 s" 1 sr" 1 (66) 
Jcb c \m x c l J 

where iV( 3 He, tx) and N(D } Tx) are given in Figs. 7(a) and 7(b). 

Having only one free parameter (baryonic density), the standard big bang nucleosyn- 
thesis is in beautiful agreement with the observations of 4 He, 3 He, D and Li, as well as 
with 3 neutrino flavors observed at LEP. Hence, our cascade nucleosynthesis calculations 
give upper limits to any hypothetical high energy process at 10 3 < z < 10 6 . However, 
one should not forget about the potential for cascade nucleosynthesis to give rise to small 
corrections for standard big bang nucleosynthesis in the production of 3 He and D. In order 
to assist such small corrections, we show in Fig. 9 the ratio of D to 3 He production by 
cascade nucleosynthesis as a function of both z c and Tx ■ 
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Figure 1: Mean interaction rates at for all of the cascade processes considered in the 
present work. Interactions with radiation at z = are shown by heavy lines: photon- 
photon pair production (full line); inverse Compton scattering (dot-dash line); photon- 
photon scattering (dotted line). For other epochs, the interaction rates for particles of 
energy E are obtained by reading the rates at energy (1 + z)E, i.e. by shifting the corre- 
sponding curve by a factor (1 -\- z) to lower energies. Interactions with fully ionized matter 
are shown by the thin lines: ordinary pair production (long-dashed line); bremsstrahlung 
(full line); Compton scattering (short-dashed line); pion photoproduction (dot-dot-dot- 
dash line). Note: we assume h = 0.7 and = 0.025 for interactions in matter; for this 
plot we use fully ionized matter (even though matter is neutral at z = 0) because, at most 
redshifts we consider in this paper, the matter is fully ionized. 



20 



S 



O 




-2-10 1 
log(£/GeV) 



o 

I 

> 

CD 

o 




S 20 



-2-10 1 
log(£/GeV) 



Figure 2: (a) The effective cross sections used for photodisintegration of 4 He into 3 He 
(thin histogram) and D (thick histogram). We estimate that the effective cross sections 
for disintegration of 4 He during photoproduction into 3 He and into D are approximately 
equal; the effective cross section for either process is shown by the dotted histogram 
labelled "7r". (b) The cross sections of part (a) multiplied by E~ ' 7 ; in the cases of 
photodisintegration into D, and disintegration of 4 He during photoproduction, the curves 
have been multiplied by 10. 
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Figure 3: The pair production cross sections for hydrogen (lower curves) and helium 
(upper curves) for the two cases: fully ionized matter (full curves), and neutral matter 
(dotted curves). 
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Figure 4: The escaping spectrum of photons from cascades on black body radaition for 
E m = AE C (see text) in which photon-photon scattering is included (solid curves) or 
neglected (dashed curves). The lower (thick) curves show results from the present Monte 
Carlo/Matrix calculation while the upper (thin) curves are from Svensson and Zdziarski 
[19]. Results of Svensson and Zdziarski are arbitrarily normalized such that the case for 
no scattering has rj(E m ) = 1. Also shown are our approximate spectrum given in Eq. 17 
(dotted curve), and the spectrum used in ref. [1] (chain curve). (All spectra cut off at 
E = E m .) 
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Figure 5: The number of (a) 3 He nuclei and (b) D nuclei produced per GeV of cascade 
energy at redshift z c . Results are shown for various h (and dotted curves - 0.4 

(0.125); full curves - 0.7 (0.025); and dashed curves - 1.0 (0.01). 
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Figure 6: The number of (a) 3 He nuclei and (b) D nuclei produced per GeV of cascade 
energy at redshift z c . Results are shown for h = 0.7 and = 0.025 and different as- 
sumptions: accurate calculation including all effects; photon-photon scat- 
tering neglected; — — disintegration of 4 He during pion photoproduction neglected; 

asymptotic pair production and bremsstrahlung cross sections used, and 

disintegration of 4 He during pion photoproduction neglected; • • • approximate 

treatment. 
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Figure 7: The number of (a) 3 He nuclei and (b) D nuclei produced per GeV of cascade 
energy as a function of mean decay time. Results are shown for various h (and dotted 
curves - 0.4 (0.125); full curves - 0.7 (0.025); and dashed curves - 1.0 (0.01). 
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Figure 8: Upper limit, 0^, to the fraction of the present closure density that massive 
particles would contribute if they had not decayed with with mean decay time tx, multi- 
plied by f c b c /ftb. Results are shown for various h (and dotted curves - 0.4 (0.125); 
full curves - 0.7 (0.025); and dashed curves - 1.0 (0.01). Also shown are results we 
would obtain if we used asymptotic pair production and bremsstrahlung cross sections 
(—• — • — • — • — ), and results of Ellis et al. [19] ( — • • • — • • • — ). 
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Figure 9: Ratio of D to 3 He for massive particle decay as a function of: (a) redshift at 
decay, z c ; (b) mean decay time, tx- Results are shown for various h (and dotted 
curves - 0.4 (0.125); full curves - 0.7 (0.025); and dashed curves -1.0 (0.01). The dot-dash 
curve shows the result for h = 0.7 and = 0.025 that would be obtained if disintegration 
during photoproduction were neglected. 
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